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ABSTRACT 

Large Hll regions, with angular dimensions exceeding 10 pc, usually enclose numerous 
massive 0-stars. Stellar winds from such stars are expected to play a sizeable role in 
the dynamical, morphological and chemical evolution of the targeted nebula. Kine- 
matically, stellar winds remain hardly observable i.e., the typical expansion velocities 
of wind-blown bubbles being often confused with other dynamical processes also regu- 
larly found Hii regions. However, supersonic shock waves, developed by stellar winds, 
should favor shock excitation and leave a well-defined spectral signature in the ionized 
nebular content. In this work, the presence of stellar winds, observed through shock 
excitation, is investigated in the brightest portions of the Galactic IC 1805 nebula, a 
giant Hii region encompassing at least 10 0-stars from main-sequence 09 to giant 
and supergiant 04. The use of the imaging Fourier transform spectrometer SpIOMM 
enabled the simultaneous acquisition of the spectral information associated to the 
HaA6563 A, [Nii]AA6548, 6584 A, and [Sii]AA6716, 6731 A ionic lines. Diagnostic 
diagrams, first introduced by Sabbadin and collaborators, were used to circumscribe 
portions of the nebula likely subject to shock excitation from other areas dominated 
by photoionization. The gas compression, expected from supersonic shocks, is investi- 
gated by comparing the pre- and post-shocked material's densities computed from the 
[s !i] A6731 ^^^^ ratio. The typical ||| A6548 ^^^'^ ratio slightly exceeds the theoretical 
value of 3 expected in low-density regimes. To explain such behavior, a scenario based 
on coUisional de-excitations affecting the [N ii] A6548 A line is proposed. 
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1 INTRODUCTION 

The presence of ionized material in the interstellar medium 
(hereafter, ISM) can be attributed to two distinctive mecha- 
nisms. First, photoionization of the surrounding neutral gas 
by the strong, energetic ultra-violet (UV) flux of nearby mas- 
sive stars is largely responsible for the detection of the stan- 
dard hydrogen Balmer series, typically used for the morpho- 
logical and kinematical description of Hll regions. Secondly, 
stellar winds with high terminal velocities and violent super- 
nova blasts are commonly associated with the propagation 
of transonic and supersonic shock waves in the surrounding 
medium. The important increase of the post-shocked gas' 
temperature favors its io nization through s hock excitation. 

Pioneering work by ISabbadin et al\ |l973) has com- 
pared specific fiux ratios for a variety of ionic transitions 
in the ISM optical gas. This allowed the authors to ap- 
proximately separate standard Hll regions and planetary 
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nebulae (hereafter, PNe) mostly governed by photoioniza- 
tion from shock-dominated supernova remnants (hereafter, 
SNRs). These diagnostic diagrams were often used, in lit- 
erature, to classify large amount of ionized targets in large- 
scale objects, for example more-or-less dist ant gala^des (e.g., 
iMagrini et aZ]|2003l : [Riesgo fc Loped bood '). This has led to 
emission-line ratio plots in which each object is usually sta- 
tistically represented by a single point. Obviously, intrinsic 
variations, within a given object, can be i nvestigated by tar- 
getin g Galactic objects individu ally (e.g.. |Phillips fc Cuestal 
199d. ll999l : IPhillips et fflEll2010l l . This allows to spatially re- 
solve much smaller (^0.1 pc) structures and artifacts char- 
acterized by peculiar line ratios that are, otherwise, unre- 
vealed or statistically negligible in observations using poorer 
angular resolutions. The investigation of close ionized ob- 
jects has already revealed that photoionization and shock 
excitation can both be found in individual regions. 

Kinematically speaking, standard Hll regions encom- 
passing massive stars earlier than 06 have shown little indi- 
cation of a strong impact on the surrounding gas attributed 
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to stellar winds, referred earlier as a potential source for 
shock excitation in nebulae. This remains usually true even 
when using observations with high spectral resolutions that 
allow to measure velocity fluctuations down to a few km 
s~^. The main reason for this resides in the double-shock 
model said to accurately describe the dynam ical evolution 
of wind-blown bubbles l|Weaver et al\ Il977l ) . The reverse 
shock quickly (i.e., often within less than a spatial element 
of resolution) converts high-velocity (> 1000 km s~^) ter- 
minal winds into low-velocity, high-temperature gas. The 
thermal energy of this hot, pressurized material then ini- 
tiates and fuels the expansion of the dense shell of post- 
shock ISM material found at much greater distances with 
respect to the central star. The typical expansion veloc- 
ity of the shocked shell is usually a few km s~^ above 10 
km s~^, roughly the speed of sound for warm (~10 000 
K) gas. Unfortunately, in complex tridimensional ge- 
ometries, this kind of velocities could be easily confused 
with standar d accelerated outflows in Hll regions, such as 
Champagne jTenorio-Tagle 19791: Bodenheimer et al\ 19791) 
or ph otoevaporated ( Arthur fc Hoare ^ '2006*: 'Mellema et al! 
20061) flows, turbulent motions (.J oncas & Roy 1984; 
Arsenault fc Rovlll988l : iLagrois fc Joncasii201ll ) or fluid in- 
stabilities. This work mainly explores the possibility of de- 
tecting shock waves associated, in particular, to stellar winds 
using the emissive properties of the ionized gas rather than 
a more typical, and not always successful, approach based 
on the information retrieved from radial velocities and non- 
thermal line widths. 

The IC 1805 nebula is located in the Perseus arm of 
our Galaxy. The most massive stars of the Melotte 15 star 
cluster are currently responsible for the energetic support of 
the large Hll region. The south-central portion of the nebula 
(see Figure 1), in the direct vicinity of the st ar cluster, is gas- 
rich a nd will be used in this work to fill the lSabbadin etali 
l| 19771 ) diagrams. Our goals are to (1) obtain reliable, high 
signal-to-noise ratios (hereafter, S/N), spectra of the opti- 
cal emission in the brightest portions of the IC 1805 star- 
forming complex, (2) provide the corresponding series of 
line-ratio diagnostic diagrams, and (3) investigate the im- 
pact of photoionization and shock excitation in the targeted 
gas volume. 

We present, in § 2, the Calactic Hll region IC 1805 and 
its associated star cluster, Melotte 15. Information related to 
our spectrometric observations and methods used for data 
reduction are detailed in § 3. Results of our study and di- 
agnostic diagrams are provided in § 4. Interpretation and 
discussion follow in § 5. Summarized results and general con- 
clusions will finally be provided in § 6. 



2 THE IC 1805 GAS COMPLEX 



Introduced by IWesterhoutI l| 19581 ) as an irregular nebula of 
medium brightness, the Galactic Hll region IC 1805 is often 
referred as the Heart nebula due to the heart-shaped mor- 
phology of its southern hemisphere (60° 28' ^ ^2000 =^ 62° 10'). 
It is suggested that outflows from the OB association 
Melotte 15 (02000 = 02*^32"' 45", feooo = 61°26'40") are re- 
sponsible for the actual expansion of the large H 11 re- 
gion/superbubble l|Basu et aZ.lll999l ). Proper motion inves- 
tigation has allowed the identification of 126 stars intrinsic 



to the star cluster (|Shi fc Hull 19991 ). The authors confirmed, 
through spectroscopic observations, the presence of numer- 
ous massive stars; about forty from spectral types 04 to 
B2 (including 10 O-type stars) have been identified with 
a high probability of membership. An important propor- 
tion of these massive stars are still found on the main- 
sequence branch i ndicating a very young cluster age eval- 
uated at 2.5 Myr l|Llorente de Andres et al! ' 1982). The dy- 
namic al age of the IC 1805 ne bula being estimated at 14 
Myr l|Lagrois fc Jonca3 l2009bl ) , it appears that the large 
Hll region was formed and shaped by the mechanical de- 
posit, in the ISM, of energy (i.e., photon flux, stellar winds, 
supernovae) attributed to a succession of different star clus- 
ters, Melotte 15 being the most recently formed. No non- 
thermal emission being actually detected in the large star- 
forming complex, this suggests that the impact of old super- 
novae could be negligible in the nebula's current dynamics. 
UBV photometry, properly corrected for absorption, puts 
the Melotte 15 complex at an heliocentric distance of 2.35 

kpc (iMassev et a;.ll 1995 1. 

iLagrois fc Joncas! (|2009al lbl) provide a high spectral- 
resolution (< 15 km s~^) kinematical investigation of the 
H""" content on an appreciable fraction of the IC 1805 ionized 
extent. The northernmost portions (63°20' < (520oo 66°50') 
of the large ovoid bubble have shown kinematical evidences 
for high-latitude gas venting and the early development 
of a Galactic chimney. It appears that leaking UV pho- 
tons, emanating from the Melotte 15 star cluster and un- 
absorbed by the denser low-declination material, may par- 
ticipate in the sustainment of the high-temperature gas 
found in the Galactic halo. In the southern half of the neb- 
ula, schematically represented in Figure 1, the last resid- 
ual fragments of an old giant molecular cloud appear to 
be eroded by the UV flux and stellar winds of the nearby 
massive stars. As a result, series of criss-crossing acceler- 
ated flows a re detected in agreement with the Champagne 
phase (e.g., Tenorio-Tagle 19791: Bodenheimer et al. Jl979l : 



iTenorio-Tagle fc Bediiniri98ll ? lTenorio-Tagle et aLlll982l ). 

Ionization fronts, stellar winds, and Champagne shocks 
are therefore expected to have developed in the vicinity of 
the Melotte 15 star cluster. This could suggest that shock 
excitation may be responsible for the presence of a certain 
fraction of the observed ionized material although a well- 
defined kinematical signature attributed to shock waves has 
not been formally identified in our previous study. This will 
be largely addressed in the following sections. 



3 OBSERVATIONS AND DATA REDUCTION 

Observations of the IC 1805 optical gas complex were per- 
formed during the nights of 2008 September 24-25, 25-26 
using the Ritchey- Chretien 1.6m telescope of the Obser- 
vatoire du Mont-Megantic (OMM). The data were gath- 
ered using the imaging Fourier transform spectrometer SpT 
OMM (Spectrometre Imageur de I'Observatoire du Mont 
Meganti c). Technical details regarding SpIOMM are pro- 
vided in lBernier et al\ (l2008l ). lDrissen ef al\ l|2008h and ref- 
erences therein. 

For each pixel of the detector, the instrument has the 
capacity to obtain the emission spectrum of the correspond- 
ing nebular-gas column in selected bandwidths of the optical 
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Figure 1. Schematic diagram of ttie souttiern portions of the IC 1805 nebu la (circumscribed by the dotted line). The regular and thick 
lines respectively traces out the H+ and CO (2-1) structures (see Figure 1 of lLagrois fc Joncaj ll2009al )). The thin-line box approximately 
delineates the gas-rich region in the vicinity of the Melotte 15 star cluster. The thick-line box approximately indicates the extended FOV 
of our spectrometric observations. Also included to the diagram, the nearby H II region W3 and the supernova remnant HB3 (long-dashed 
lines). 



Table 1. Observational Parameters at Data Acquisition 



Field 



Observational Parameters 


East 


West 


Optical center 


O2000 
(52000 


02'^33™55= 
61°27'30" 


02''32™35= 
61°27'30" 


Seeing at acquisition 


(") 


2.5 


2.8 


Spatial resolution 


(" pixcl-i) 
(pc pixel ~-'^) 


2.2 
0.026 


2.2 
0.026 


Spectral Range 


(A-A) 


6480-6820 


6480-6820 


Spectral resolution 


(A) 
(km s-i) 


3 
150 


3 
150 


Exposure time per channel 


(s) 


32 


35 


Number of channels 




249 


249 



regime between 3500 A and 9000 A. Prior to entering the 
observation mode, the spectral resolution is fixed by the ob- 
server between R = 1 and 25 000. 

For this work, the use of an interference filter enabled 
the acquisition of the spectral information in the red portion 
of the electromagnetic spectrum between 6480 A and 6820 
A. This allowed to obtain, simultaneously, emission spec- 
tra displaying the HaA6563 A, [Nll]AA6548, 6584 A, and 
[Sii] AA6716, 6731 A ionic transitions. A spectral resolution 
of R = 2100 was judged appropriate corresponding to a res- 
olution of roughly 3 A or 150 km s~^ centered on the Hq 



rest frequency. The field-of-view (FOV) is roughly 12' x 12' 
while data were spatially binned 2x2 during acquisition to 
reduce readout times. 

The resulting raw data cube is formed, for each pixel, of 
a discrete interferogram. Following classical data reduction 
for CCD observations (i.e., bias subtraction and flatfield cor- 
rection), all panchromatic images were realigned to correct 
for guiding errors. According to the photometry of stars lo- 
cated in gas-depleted zones, all images were also corrected 
for variations of the sky transparency occurring during data 
acquisition. Each interferogram is then Fourier transformed 
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Figure 2. High-quality emission spectrum recovered from our spectrometric observations of the brightest parts of the IC 1805 gas 
complex. The raw data are displayed as the regular line. Gaussian fits are represented as the dotted line. The boxes indicates the inter- 
line channels where the noise level was estimated. All five ionic transitions, targeted in this work, are listed. The S/N varies from 24 to 
298 from [N ii] A6548 A to Ha A6563 A. 



and calibrated in velocities using a He-Ne laser (A6328 A) as 
frequency of reference. The data were again spatially binned 

2x2 during reduction in order to account for the 2 . 5 seeing 
measured at the telescope. The resulting spatial resolution 

is of 2 . 2 pixel for 325 x 335 pixel . Using the adopted dis- 
tance to Melotte 15 (see § 2) , this corresponds to an angular 
resolution of 0.026 pc pixeP^ as seen on the plane of the 
sky. 

Two cubes were obtained in the vicinity of the Melotte 
15 star cluster. The observational parameters for both cubes, 
labeled East and West, are summarized in Table 1. The he- 
liocentric correction being very similar for both cubes, the 

two were mosaicked to a new extended FOV of 21.5 x 12.5 
(see Figure 1) and over 200 000 emission spectra. Following 
subtraction of the continuum, a multi-component Gaussian 
fit procedure, written in idl, was applied to the mosaicked 
cube. The fitting procedure returns, for each component, the 
peak intensity in Analog-to-Digital Units (ADU), the line 
centroid (or velocity), and the line dispersion {a= ) . 

Since this work is most ly based on line rat ios, we followed 
the recommendation of iRola fc Pela^ ()l994h by imposing a 
S /N greater than 6 for a given Gaussian fit to be usable. For 
each spectrum, the noise level was estimated from contin- 
uum fiuctuations measured in empty channels. 

Figure 2 shows a high-quality spectrum obtained in the 
brightest parts of our mosaicked cube. This kind of spec- 



trum is typical only for the regions of high emissivity of 
our FOV. Emission-line profiles corresponding to weaker ar- 
eas are much noisier and complementary lines, found in the 
wavelength interval covered by the interference filter (see 
above) , such as the [N l] A6500 A and [He l] A6678 A transi- 
tions are not discernible from noise fluctuations. The weak- 
est line (i.e., [N ll] A6548 A) in the example provided by Fig- 
ure 2 has a S/N of 24. 

4 RESULTS 

4.1 Monochromatic maps 

Figures 3 to 7 provide the monochromatic, peak-intensity 
maps for all five emission lines investigated in this work. 
North is up and East is left in all figures. Figure 8 spec- 
tacularly displays the ^^^'^^^^gy ratio throughout the whole 
FOV (see caption). Red emphasizes regions dominated by 
Ha while bluer shades pinpoint a stronger emission of the 
[Nil] A6584 A transition (with respect to the redder areas). 
Three stars with strong Ha lines are detected (the third be- 
ing barely perceptible very close to the center of the figure) . 
At the position of all other stars, bad pixels were removed 
and replaced by mean values statistically representing the 
surrounding nebular-gas content. This figure is particularly 
useful to describe ionized features in our FOV. 
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Figure 3. Peak-intensity map of tlie [N ll] A6548 A ionic transition in IC 1805. Only emiss ion-line profiles with S/N > 6 were considered. 
A logiQ-scale filter was applied in order to emphasize fainter ionized structures. As listed bv lMassev et all lll995l ). the filled circles indicate 
the position of the six O-stars found in the FOV. Dotted contours are traced out using an integrated intensity map of the molecular 
emission obtained from the FCRAO CO (1-0) observations. The spectral collapse of the molecular data cube was made for local standard 
of rest (LSR) radial velocities bet ween -22.90 and -63.30 km , approximately corresponding to the extrema in velocities measured in 
the H+ gas closest to Melotte 15 jLaerois &: Jonca3l2009al ). Eleven contour levels are drawn between 2 and 12 K km s'^ The FCRAO 
CO (1-0) subcube has a noise level slightly below 1 K. 




Figure 4. Peak-intensity map of the Ho A6563 A ionic transition in IC 1805. More details are provided in the caption of Figure 3. 
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Figure 5. Peak-intensity map of the [Nil] A6584 A ionic transition in IC 1805. More details are provided in the caption of Figure 3. 
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Figure 6. Peak-intensity map of the [Sll] A6716 A ionic transition in IC 1805. More details are provided in the caption of Figure 3. 



The Five College Radio Astronomy Observatory 
FCRAO) CO(l-O ) survey of the 2'"* Galactic quadrant 
Hever et fflZ.lll99§ ) has revealed very faint emission corre- 
sponding to tenuous molecular material at the position of 
the bright, central ionized structure in Figures 3 to 7. This 
suggests the presence of a large molecular fragment, sur- 
rounded by the most massive stars of the Melotte 15 clus- 
ter, that has undergone almost full erosion by the UV flux 



and stellar winds of the nearby ionizing sources. Figure 8 re- 
veals the filamentary nature of the central structure. These 
filaments most likely trace out the spatial disposition, on 
the plane of the sky, of the eroded (sometimes fully ionized) 
molecular envelopes. The central structure is surrounded by 
diffuse ionized material likely associated to photoevaporated 
flows kinematically in agre ement with the Champagne phase 
(jLagrois fc Joncaa l2009al ). One of these flows (displaying 
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Figure 7. Peak-intensity map of the [Sll] A6731 A ionic transition in IC 1805. More details are provided in the caption of Figure 3. 




Figure 8. Image of the ;^g5g3 ratio in central IC 1805. Red shades indicate that the Ha component clearly dominates the signal. Blue 
shades show the position of particularly bright [Nil] A6584 A emission. Stars were removed if not for three of them having well-defined, 
strong Ha emission (see text). This image was not processed using the information retrieved from the Gaussian fits but rather by simply 
collapsing the initial mosaicked cube in the channel intervals respectively enclosing the Ha A6563 A and [Nil] A6584 A ionic lines (e.g., 
see Figure 2). 



blue shades in Figure 8) is particularly well-defined, prop- 
agating from the bottom-center of our FOV approximately 
toward the upper-left corner of Figure 8. West of the cen- 
tral structure, the spatial arrangement of the nebular gases 
is extremely complex. Ionized filaments are still perceptible 
although the contrast with the diffuse surroundings is not 



as clear. Toward the western boundary of Figure 8, even the 
Hq emission appears darker, obscured by interstellar dust. 
Shock excitation in the vicinity of the bright, central ion- 
ized structure, displayed in Figures 3 to 8, will be largely 
discussed and investigated in § 5.2.1. 

Toward the south-eastern portion of our FOV, a small 
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Figure 9. Image of the " [ n nfAessT**^^ ratio in central IC 1805. Blue shades indicate that the [Nil] component clearly dominates the 
signal. Red shades show the position of particularly bright [S ll] emission. The image was processed identically to Figure 8. 



CO fragment is clearly detected in the FCRAO survey at 
the systemic velocity of the IC 1805 region (see CO con- 
tours in Figures 3 to 7). Its ionized counterpart is revealed 
by a thin, rounded ionization front especially visible in the 
lower-left corner of Figure 8 (see also Figure 19a). This spe- 
cific shape of the ionization front results from the vast ma- 
jority of the ionizing sources in Melotte 15 being located 
behind th e molecular clump. This was kinematically con- 
firmed, in lLagrois fc Joncad (|2009al ). by the detection of an 
accelerated ionized outfiow moving away from the observer 
(see § 5.2.2.1). This fiow eventually collides or simply coin- 
cides in lines-of-sight with the south-central/north-east fiow 
mentioned in the previous paragraph. This results in partic- 
ularly complicated kinematical motions in central IC 1805 
with Ha non-thermal line widths approaching the supersonic 
regim e (^10 km s~^) accordin g to high-resolution observa- 
tions (|Lagrois fc Jonca3l2009al ). 

In the same area of Figure 8 (see also Figure 19c), a 
cylindrical, cigar-like feature appears to be associated with 
an isolated star with strong Ha lines (see above). The fea- 
ture has bright H+ and rims but is almost completely 
gas-deprived near its center. This structure as well as the 
rounded ionization front found in the lower-left corner of 
Figure 8 will be discussed in § 5.2.2. 

Figure 9 was processed identically to Figure 8 and dis- 
plays the "[NnflesM'^'^^ ratio (where the numerator is the 
sum of both lines of the [Sll] doublet). Again, blue shades 
identify nitrogen-rich zones while red areas indicate that a 
sizeable role is played by sulfur in the overall gas emissiv- 
ity. One immediately notices the complexity of the chemical 
properties in IC 1805. While [N ll] clearly dominates over 
[Sll] where Ha is the brightest, weaker zones (in Ha in- 
tegrated intensity) show an important increase of the rel- 
ative contribution of the sulfur material. In particular, the 
reader's attention is directed toward the north-eastern fil- 



ament of the bright, central structure described above. Se- 
ries of small, quasi-circular blobs are found. The top three 
reveal two nitrogen-dominated features (in blue) and one 
with a stronger relative contribution in [Sll] (in red). In- 
vestigating their corresponding spectrum, the [S ll] intensity 
remains roughly constant from one blob to another while 
the [Nil] lines sufTer a strong attenuation along the line-of- 
sight of the "redder" one (which explains the relatively poor 
contribution of nitrogen towards it). This behavior is not 
generalized to our whole FOV although similar features can 
be found here and there. 



The supernova remnant HB3 (see Figure 1) is located 
too far away from our FOV to have had a sizeable im- 
pact on the chemical properties in central IC 1805. Alter- 
natively, we can argue that the old, large molecular cloud, 
that gave birth to Melotte 15 roughly 2.5 Myr ago, may 
have had an inhomogeneous distribution of its chemical 
compounds. Old supernovae in IC 1805, whose non-thermal 
emission has vanished since, could also be held responsi- 
ble. We reiterate that the IC 1805 region was most likely 
formed by a succession of different star clusters (see § 2) 
i.e., although no indication for supernova remnants is cur- 
rently found inside the large Hll region, it is highly prob- 
able that the inner zones of IC 1805 were, at some time in 
the past, disturbed by supernova events associated to pre- 
vious generations of massive stars. These stars could have 
had an intrinsic inhomogeneity in their inner nitrogen dis- 
tribution which eventually led to an anisotropic dispersion 
of these chemical compounds, products of the CNO cycle, 
as each s tellar object reach e d the e nd of its life (e.g., see the 
works bv lMac Alpine et al\ ||2007|). MacAlpine fc Satterfieldl 
(' 20081 ') ■ and lCharlebois et al\ l|2010l ') on the Crab nebula). 
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4.2 Diagnostic line ratios 

4-2.1 Definition 

Line ratios, in the literature, usually correspond to ratios 
between two (or more) line fluxes. For a given emission line, 
the line flux is proportional to the product between its peak 
intensity and its width as returned by the Gaussian flt. In 
this work, the relatively low spectral resolution used (see 
§ 3) is roughly a factor 10 to 20 greater than typical non- 
thermal velocity fluctuations, along the line-of-sight, found 
in IC 1805 using high-resolution observations of the kine- 
matics ([Lagrois & Joncas 2009a). Hence, in this work, the 
width of each emission line is entirely dominated by the in- 
strumental response and the returned widths are very sim- 
ilar, from one ion to another, independently of the position 
in the FOV. Therefore, line ratios, in the following discus- 
sion, were strictly estimated using returned values for peak 
intensities while line widths were simply not considered. 

4-2.2 Line profile selection 

To assure the selection of the most reliable emission-line 
proflles in our sample of over 200 000 spectra, a series of 
conditions are proposed. All of them need to be fulfilled 
otherwise the given spectrum is rejected from the discussion 
to follow. The conditions are summarized as follow: 

(i) The profile must show all five lines (i.e., HqA6563 A, 
[N ii] AA6548, 6584 A, and [S ii] AA6716, 6731 A) investigated 
in this work with reliable S/N above 6. 

Two conditions are proposed in order to confirm that 
the kinematical properties of a given profile are physi- 
cally acceptabl43- More specifically, the [N ii] A6548 A and 
[N ii] A6584 A transitions must exhibit identical non-thermal 
motions since both lines are emitted by the same ions. 
The same argument also applies to the [S ll] A6716 A and 
[S ii] A6731 A lines. 

(ii) First, the difference, in centroid velocities, between 
both lines of a given doublet (i.e., [Nil] or [Sll]) must not 
exceed 15 km s~^. This corresponds roughly to 1/lOth the 
spectral resolution of our observations (see § 3) and there- 
fore roughly accounts for the uncertainties on the Gaussian 
fits. Centroid- velocity differences, along a given line-of-sight, 
between the H^, N"*", and S"*" material are allowed (and ex- 
pected!) since all three ions are not necessarily co-spatial in 
central IG 1805. 

(iii) Secondly, the difference, in measured line widths, be- 
tween all five components must not exceed 15 km s~^. This is 
expected from all lines being dominated by the instrumental 
response (see § 4.2.1). If observations with high spectral res- 
olution (e.g., of the order of a few km s~^) were used here, a 
condition similar to the previous one (on centroid velocities) 

^ Although this paper is based on line-ratio measurements, kine- 
matical properties were also carefully investigated. The ques- 
tion arises: could the line fluxes of kinematically "unaccept- 
able" emission- line profiles (i.e., that do not fulfill Conditions 
(ii) and/or (iii)) be trusted? At this point, we have no reliable 
method to further investigate the question. We therefore chose 
to avoid any possible controversy /confusion on the matter and 
simply reject those kinematically peculiar spectra. 



would have been required (i.e., identical line widths, within 
the uncertainty bars of the Gaussian fits, between both lines 
of a given doublet, [Nil] or [Sll]). 

Finally, specific line ratios of the [N ll] and [S ll] transi- 
tions must agree with standard theoretical models of ionized 
nebulae. 

(iv) The computation of the [|"| a673i ^^'^^ ratio must lead 
to a flnite electron density. The upper and lower uncertain- 
ties on the density measurement (see § 4.2.3) must also be 
finite. 

(v) The computation of the || j A6548 ^^'^^ ratio must not 
largely deviate from the theoretical value approaching 3 in 
low-density regimes typically found for Hll regions. Small 
deviations from 3 are expected (see § 4.2.4) and explained 
(see § 5.1) although values ^2 or 3>4 would seem particu- 
larly hard to reconcile with the theory and would therefore 
demand to be rejected. 

All conditions taken in consideration, only 3 057 
emission-line profiles were retained, most of them being as- 
sociated to the bright, central structure found in the over- 
lapping region between the eastern and western field (see 
§ 5.2.1). Obviously, the summation of both cubes largely 
contributes to increase the data quality for duplicated pix- 
els not only by increasing the peak signal for all lines but 
also by flattening noise fluctuations in empty channels (see 
Figure 2). 

Considering the very large number of emission-line pro- 
files available in our initial data set (more than 200 000), 
we are fully aware that the conditions listed above drasti- 
cally reduce the size of the retained sample. On the other 
hand, these conditions assure that the retained profiles are 
undoubtedly the most reliable available. 

We reiterate that a S/N greater than 6 is required to 
appropriately conduct an investigation on line ratios. How- 
ever, a "physic al detection" [of a given line] is considered 
when S/N>3 (|Rola fc Pelatlll994l '). This wiU be used in a 
later subsection to investigate notable, although not neces- 
sarily reliable, line ratios for particular structures in IC 1805 
(see § 5.2.2.2). 

4-2-3 Electron densities computed from the [S u] doublet 

The [g "j A6731 ^''^^ ratios were computed and used as input 
value s for the Fivel FORTRAN procedure (|De Robertis et aZI 
Il987h adapted, for convenience, in idl. Assuming a con- 
stant electron temperatu re of 7 400 K throughout IC 1805 
(|Lagrois fc Joncaslf2009ah . mean electron densities (p), in 
the S"*" ionic volume, were obtained for each nebular-gas col- 
umn (i.e., line-of-sight). 

Since uncertainties cannot be recovered from the use of 
the Fivel procedure, error bars on density measurements p 
were computed as follows. For each line of the [Sll] dou- 
blet, the uncertainty on the peak intensity, as returned by 
the Gaussian fi t proce dure, is provided by Equation 4a of 
[Landman et al\ (|l982l '). From there, the statistical error on 
each line ratio is simply provided by the standard propa- 
gation of uncertainty. For each ratio, minimal and maximal 
plausible values are then easily obtained by respectively sub- 
tracting and adding this calculated uncertainty. Using these 
lower and upper limits on ||||] ^fyg^ , maximal (pmax) and 
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Figure 10. Panel (a): Histograms of the mean electron densities (p) in IC 1805 computed from tlie | g "j line ratio. The black-line 

distribution corresponds to the 3 057 emission-line profiles fulfilling all five conditions of § 4.2.2. The blue curve is a subset of the black 
distribution, isolating the electron densities of specific points suggesting shock excitation. The red distribution uses the same points as 
the blue histogram although densities were computed from the post-subtraction [S ll] doublets. The two colored histograms are discussed 
in § 5.2.1 (see text). Panel (b): Spatial distribution of all 3 057 electron- density measurements forming the black histogram of Panel (a). 
Densities between 1 and 75 cm~^ (purple dots), 75 and 150 cm"'' (blue dots), 150 and 200 cm"'' (green dots), 200 and 300 cm"'' (orange 
dots), and 3004- cm^'' (red dots) are displayed. 



minimal (pmin) densities, associated to p, can be obtained 
via Fivel. An asymmetrical error bar usually results where 
the upper (lower) uncertainty is provided by the difference 
between pmax (p) and p {pmin)- Only p values where both 
Pmin and pmax are finite were retained (see § 4.2.2). Ac- 
cording to Fivel, this implies that pmin and pmax must be 
greater than 10 cm~^ and lower than 15 000 cm""' respec- 
tively (these two values mostly depend on our choice of a 
constant electron temperature of 7400 K in IC1805). Our 
method considers that the statistical error on p is entirely 
dominated by the data quality (i.e., S/N of the [Sll] lines) 
and that the uncertainties on the atomic parameters, used 
by Fivel, are negligible. However, using pmax and pmin to cal- 
culate the statistical error on p probably overestimates the 
actual uncertainty on the electron density measurement. 

The distribution of the retained electron densities is 
provided by the black curve in Figure 10 a (the blue and red 
distributions will be discussed in § 5.2.1). The histogram has 
a mean of 155^ 75^ cm~^ where the error bars correspond to 
the mean upper ((pmax — p)) and lower ((p — pmin)) uncer- 
tainty retrieved from the 3 057 electron-density points pre- 
served. Panel (b) provides the spatial distribution of these 
points, superimposed to a black-and-white reproduction of 



Figure 4. Purple dots have the lowest density values while 
red dots, the largest (see caption). 

4.24 The [Nil] line ratio 

Figure 11 provides the distribution of the ||j ^g54g hne ra- 
tio in the vicinity of the Melotte 15 star cluster. The vertical 
long-dashed line indicates the theoretical value as suggested 
by the ratio of the transition probabilities of both ionic 
lines in a low-density regime l|Osterbrock fc Ferlandl |2006| . 
Chapter 5). All 3 057 points mentioned in § 4.2.3 formed 
the regular-line distribution. As shown by the thick-line his- 
togram, the presence of large values i.e., above ^3.5—3.6, is 
reduced by considering only data points of very high quality 
(i.e., requiring here, arbitrarily, that S/N > 20 in the first 
condition of § 4.2.2). Uncertainties on the mean and one- 
standard deviation of both histograms are statis tical uncer- 
tainties retrieved assuming normal distributions (|Bevingtonl 
1969, Chapter 5). Extreme values for ||j ^gg4g i.e., below 
2 or well-above 4, were not measured. Hence, Condition 5, 
listed in § 4.2.2, was not used in this work to reject partic- 
ular emission-line profiles although it could used by other 
authors for future, similar studies. 
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Figure 11. Histograms of the ^6543 ratios in IC 1805. The theoretical value of 2.94, obtained in a low-density regime, is 

indicated by the long-dashed line. The regular-line distribution corresponds to the 3 057 emission-line profiles fulfilling all five conditions 
of § 4.2.2. The thick-line distribution considers only emission-line profiles with S/N greater than 20 (see text). Statistical properties of 
both histograms are provided within the figure itself. 



Considering the relatively high S /N used here (above 20 
for the thick-line histogram of Figure 11 hence suggesting 
a very good signal quality), deviations from ^g54g ~ 3 
seem genuine and will be addressed in § 5.1. 



4-2.5 Line-ratio diagnostic diagrams 

Panels (a) of Figures 12 to 14 provide a series of line-ratio 
diagnostic diagrams commonly used in optical observations. 
Each diagram contains the 3 057 points already discussed 
in §§ 4.2.3 and 4.2.4. Areas labeled "Hll regions", "SNRs", 
a nd "PNe" vyere a l l repr oduced according to Figures 1 to 3 
of lSabbadin et al\ (|l97'7 ). Long-dashed lines, in Figures 13tt 
and 14a indicate the lower and upper limits for [g"| a673i 
that provide finite values of electron densities p given a con- 
stant temperature of 7400 K in IC 1805 (see § 4.2.3). 

In all three diagrams, crosses, colored in blue, are used 
to represent points found (in majority) within the "Hll re- 
gions" area. The relatively high numbers of crosses is not 
surprising for IC 1805, already cataloged as an Hll region. 
A certain fraction of the sample, however, falls outside this 
area. These points are symbolized as red filled circles and 
point approximately toward the "SNRs" area, hence sug- 
gesting evidence for shock excitation in the targeted ISM 
volume. This will be discussed in § 5.2. Out of the 3 057 
points retained for this study, respectively 134 (4.4%), 67 
(2.2%), and 84 (2.7%) are displayed as red filled circles in 
Figures 12a to 14a. 

Panels (b), in all three figures, provide the spatial distri- 
bution of all crosses (shown as blue dots) and circles (shown 



as red dots) displayed in the diagnostic diagram of their re- 
spective Panel (a). As in Figure 106, a black-and-white re- 
production of Figure 4 was used. Red dots in Figures 126 to 
146 all point at the same areas and reveal that the southern 
portions of the bright, central structure (^2000 <61°26'30") 
may encompassed material subject to ionization by shocks. 



5 DISCUSSION 

5.1 The [Nil] line ratio: 
theoretical value 



deviations from the 



As demonstrated in § 4.2.4, values for the [^"j^"^^^ fine 

' [N II] A6548 

ratio above 3.5 and approaching 4 are partially attributed 
to poorer data quality especially affecting the lower-signal, 
noisy [Nil] A6548 A line (e.g., see Figure 2). Figure 11 indi- 
cates that the use of profiles whose emission lines are strictly 
characterized by "very" high S/N (in our case, greater than 
20) contributes to attenuate the right tail of the distribu- 
tion. The asymmetry coefficient (i.e., skewness) is reduced 
from roughly 0.4 to 0.2. Consequently, the width of the dis- 
tribution is statistically narrower while the mean of the his- 
togram slightly varied toward the expected value of ~3 for 

[Nil] A6584 
[Nil] A6548 ■ 

Nonetheless, independently of the threshold in S /N used 
(e.g., 6 or 20), Figure 11 reveals that the mean []^||| xtslt 
ratio exceeds, by a few tenths, the theoretical value. We 
propose that this could be attributed to a simple statistical 
effect by which the ratio of transition probabilities between 
the [N 11] A6548 A and [N 11] A6584 A transitions does not 
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Figure 12. Panel (a): Diagnostic diagram of the logio /(^ [ s n ] ^) ] [ i(|n n] ) ] ''elation. By definition, I(Ho) is the peak intensity 

of the Ha A6563 A line while I([Sll]) and I([Nll]) correspond respectively to the sum of both lines' peak intensity of the [Sll] and [Nil] 
doublet. While blue crosses indicate photoionization effects typically expected in Hll regions, red filled circles deviate toward the shock- 
dominated regime of SNRs (see text). Panel (b): Spatial distribution of all 3 057 points found in Panel (a). Blue crosses are represented 
by blue dots and red filled circles, by red dots. 



necessarily equal 3 given the physical properties governing 
the IC 1805 nebula. Note that the curve of transmission of 
the interference filter used for data acquisition has not re- 
vealed any significant difference between the transmission 
coefficients of both lines of the [Nil] doublet. 

Table 3.15 of lOsterbrock fc Ferlandl (|2006l ') provides, 
for different ionic transitions, the critical electron density 
above which coUisional de-excitations, caused by free elec- 
trons, become important. This density is estimated at 80 and 
310 cm"^ respectively for the [Nil] A6548 A and [Nil] A6584 
A transitions. Hence, for a given electron density between 
these two values and characterizing the N"*" volume, the 
[N ll] A6548 A line intensity will be less than expected, a 
fraction of its flux being "lost" in (non-radiative) coUisional 
de-excitations. On the other hand, the [N ll] A6584 A line in- 
tensity will not be affected until densities reach the 310 cm~"^ 
plateau. Figure 10a provides reasonable evidences for elec- 
tron densities between 80 and 310 cm~^ in IC 1805. Hence, 
a mean ^^543 ratio of 3.3—3.4 could be explained by 
a [N 11] A6548 A line partially experiencing coUisional de- 
excitations. 

From the results presented in §§ 4.2.3 and 4.2.4 (for 
S/N>6), Figure 15 provides the electron densities (as 
measured from the ^gyg^ line ratio; see Figure 10) 



vs. the '^ce/c line ratio (see Figure 11) relation. Each 



[Nil] A6548 

_ [Nil] 
P [NlllA6548 

bin of 10 cm"'' between 10 and 500 cm 



[Nil] A6584 p^j^ represented by a black cross. For each 



the mean 

[nII] A6548 ^^^^ ratio is displayed as a red filled square. A 
standard linear regression (|Bevingtonlll969l . Chapter 4) was 
applied to the fiUed-square sample, in the 80—310 cm~^ in- 
terval, and is displayed as the dashed green line. The value of 
each bin was adequately weighted in the regression i.e., bins 
with large numbers of crosses were given a more important 
statistical weight. 

A clear, monotonic increase of the !^"! '^'^^^o line ra- 

' [N 11] A6548 

tio with increasing densities would be expected as the 
[N 11] A6548 A line becomes more-and-more affected by col- 
lisional de-excitations (i.e., decreases in intensity). As ex- 
pected, the linear relation shows a positive slope although 
its value, very close to 0, largely suggests no correlation at 
all. Indeed, the quality of the linear fit is relatively poor with 
a correlation coefficient of 0.35. Hence, we believe that no 
conclusion can be drawn from Figure 15. Two possibilities 
are however provided in order to explain such behavior. 

First, the reader should note that the uncertainties are 
relatively large for both parameters investigated. The mean 
uncertainty for the [^||] ^6548 ^^^^ ratio is 0.17 while typical 
uncertainties for the electron-density measurements p are 
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Figure 13. Panel (a): Diagnostic diagram of the logio ^ ^[g n ] ^) ] [ s A6731 ratio relation. The right-hand ordinate indicates the 
correspondence between [S ll] line ratios and electron densities p as provided by our IDL adaptation of the Fivel FORTRAN procedure. 
Panel (b): Spatial distribution of all preserved points. More details are provided in the caption of Figure 12. 



provided in the last paragraph of § 4.2.3. The poor corre- 
lation of the linear fit in Figure 15 could be attributed to 
large error bars. 

Moreover, the absence of a well-defined correlation be- 
tween p and ^gg4g would not be surprising if the N'^ 
and S"*" material are not perfectly cospatial in the vicinity of 
Melotte 15. Given a difference of more than 4 eV between 
the ionization potentials of neutral nitrogen and sulfur, it 
can be assumed that the electron densities, computed from 
the [Sll] doublet, may not refiect accurately the physical 
conditions prevailing in the ionic volume. Therefore, the 
P [Nn] A6548 diagram of Figure 15 could be an interesting 
diagnostic tool although our data set may not provide the 
appropriate spectral information to reliably construct such 
a diagram. Alternatively, given the first ionization potential 
of nitrogen (14.5 eV) and oxygen (13.6 eV), the [o||j A3729 
line ratio would likely give a more accurate indication of the 
density behavior inside the volume. 

5.2 Shock excitation in IC 1805 

Panels (a) of Figures 12 to 14 reveal that a large fraction of 
the nebular-gas content in the vicinity of the Melotte 15 star 
cluster appears to be dominated by photoionization, typical 
of standard Hll regions. However, points identified outside 
the expected regime for Hll regions and located near the 
"SNRs" area demand to be closely investigated (see § 4.2.5). 



In all three figures, these are symbolized as red filled circles 
and potentially indicate the presence of shocks in the tar- 
geted nebular volume. These areas of our FOV, suggesting 
shock excitation, are identified as red dots in all three Panels 
(b) of the same figures. 

However, at this point, shock excitation in IC 1805 re- 
mains highly hypothetical. This is somewhat surprising con- 
sidering the presence of relatively massive stars in Melotte 
15 and their associated strong stellar winds (see § 2). In or- 
der to deeply investigate the impact of shock excitation in 
IC 1805, we assume that, if present, the shock-excited ion- 
ized material is relatively confined and localized along the 
line-of-sight i.e., consi dering the imposing dime nsions of the 
large Hll region fe.g.. iLagrois fc Joncas|[2009al lbl). the inter- 
face between compressive shocks and the surrounding ISM 
could be several times smaller than the size of the nebula 
itself. This said, the expected spectral signature of shock 
excitation could be strongly diluted by photoionized fore- 
ground/background material. 

5.2.1 Central structure 

A total of 2 378 out of 3 057 emission- line profiles retained 
for this study are directly associated to the bright, central 
structure clearly visible in Figure 8. The strong emission 
of its ionized component (see peak intensities in Figures 3 
to 7) and therefore the quality of the gathered signal make 
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this structure, surrounded by the most massive stars of the 
Melotte 15 cluster, the ideal feature in our quest for shock 
excitation in central IC 1805. 

First, a series of four weaker portions [in gas emis- 
sion], specifically surrounding the central structure, were 
selected. Each of these portions was spatially binned into 
a single 1 x 1 x 249 (spatial x spatial x channels) emission- 
line profilfl For each of the 2 378 spectra selected and 
investigated in this subsection, the corresponding fore- 
ground/background emission was approximately recovered 
using a linear combination of these four 1 x 1 x 249 pro- 
files. Each of the four linear coefficients was weighted via 
the inverse of the distance separating the targeted pixel 
from the center of the corresponding weaker portion (e.g., 
the larger this distance is, the smaller is the corresponding 
linear coefficient and, therefore, the smaller is the statisti- 
cal impact of this weaker zone on the computation of the 
foreground/background spectral signature at the position of 
this given pixel). Hence, a unique foreground/background 
spectrum is constructed for each of the 2 378 points (see 
below) although neighbor pixels have [as it should be] fore- 
ground/background material with very similar spectral sig- 
natures (i.e., roughly identical linear coefficients). 

The subtraction of the foreground/background emission 

^ The spatial binning of numerous emission-line profiles into a 
single one flattens out noise fluctuations. Therefore, the resulting 
spectrum is of very high quality. 



obviously had a certain incidence on the S/N of the resulting 
profiles, the peak intensity of all five lines being reduced in 
the subtracting process. Still following all conditions listed 
in § 4.2.2, 2 229 emission-line profiles out of 2 378 were said 
to "survive" the procedure. Panels (a) and (b) of Figures 

16 present the logio [rfif^] vs. logio diagnostic 
diagram for these 2 229 points respectively before and af- 
ter the subtraction of the foreground/background material. 
Figure 16a is obviously a subset of Figure 12a and uses 
the same symbol definitions (see § 4.2.5). Note that each 
point preserves its original symbol, shape (cross/circle) and 
color(blue/red), from Figure 16a to 166. 

In Figure 16a, 2 109 are displayed as blue crosses while 
120 are symbolized as red filled circles, being located close 
but outside the lower-left corner of the "Hll regions" area. 
Note that 134 points were displayed as red circles in Figure 
12a (see § 4.2.5): 14 of those points were not considered 
here, either not associated with the central structure or have 
simply failed to fulfill the conditions of § 4.2.2 following the 
subtracting process discussed above. 

In Figure 16b, a tail precisely directed toward the 
"SNRs" area appears as numerous points migrate out of the 
"Hll regions" area. Recent surveys (e.g., iRiesgo &: Lopea 
^2006; Frew fc Parker 2 010,) have shown that the boundaries 
of the different areas circumscribed in the diagnostic dia- 
grams used here could be sli ghtly modified with re spect to 
the first work carried on bv ISabbadin et oZl (|l977l ). It ap- 
pears that the lower and upper limits, in the parameter 
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Figure 15. Diagram of the mean electron density (p) vs. ||j ^6548 ^^^^ ratio relation. The horizontal dotted line indicates the theoretical 



value of the 
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ratio retrieved 



[ N ii] A6548 predicted in a low-density regime. The horizontal regular line marks the mean A6548 

from Figure 11 (regular-line distribution). The vertical regular lines isolate the 80—310 cm~^ interval in densities where collisional de- 
excitations are expected to impact the [Nil] A6548 A line while the [Nil] A6584 A line remains theoretically unaffected (see text). Red 



filled squares correspond to the mean 



[Nil] A6584 



[Nil] A6548 '^^^'^ observed for each bin of 10 cm ^ in densities. The dashed green line indicates the 
linear fit to the sample of filled squares found between 80 and 310 cm~^. Large fluctuations of the filled-square symbols at high densities 
(/9>350 cm~'^) result from low amounts of points (black crosses) within each bin. 



space, for each regime are still uncertain to this day. We 
therefore chose to arbitrarily redefine (or expand) the zone 
of shock-excitation for each relation investigated. For ex- 
ample, in the logio [i^^f^jj] vs. logio [ /(j""]) ] relation, the 
hatched area in the lower-left corner of Figure 166 was, from 
now on, used to separate gas columns (i.e., pixels) likely 
enclosing post-shocked material from those for which pho- 
toionization effects appear to primarily dominate. Now, 509 
points out of 2 229 suggest evidence for potential shock ex- 
citation in the direct vicinity of the Melotte 15 cluster. The 
tail is formed of both blue crosses and red circles indicating 
that Figure 16 a alone would have not be sufficient in order 
to precisely pinpoint those sections of the central structure 
revealing (or hiding) a well-defined signature of ionization by 
shocks. Panel (a) of Figure 17 provides the post-subtraction 
spatial distribution of all points found in Figure 166. Points 
found in the hatched area of the lower-left portion of the 
diagram are colored in red. Others, still associated to the 
"H 11 regions" area hence suggesting photoionization effects, 
are in blue. 

Similar exercises were carried on using the correspond- 
ing subsets for Figures 13o and 14a although we will 
simply provide the results here. The zone of shock ex- 
citation was, respectively for both figures, redefined by 



l°gio [loff^D] ^ ^^'^ ^°Sio [ i(|n?i]) ] < 0.2, and also bor- 
dered by the two long-dashed lines (i.e., lower and upper 



1- ., [Sii]A6716n 
limits on [sii]A673l )- 



For the logio [i^l^] ^s. diagram, while only 

54 out of 2 229 points were being found initially outside 
the "If 11 regions" area prior to the subtraction of the 
foreground/background material, 509 points have shown 

logio [ ilisu]) ] <0-3 afterwards. Again, as in Figure 166, evi- 
dences for shocks appear for both blue crosses and red filled 
circles. Panel (b) of Figure 17 gives the spatial distribution 
of the post-subtraction points. One clearly sees the perfect 
correspondence between Panels (a) and (b) i.e., using two 
distinct diagrams the same points, that suggest shock exci- 
tation, are extracted. 

For the logio [i^jf^] }s "] a673i relation, the situ- 
ation is more complex as the "SNRs" and "Hll regions" 
areas overlap in the diagnostic tools (see Figure 14 a). The 
removal of the foreground/background material has led to 
1 419 points displaying logio [^ /(^Nnj) ] values below 0.2 al- 
though roughly half of these remain very close to the "Ifll 



relation. 



regions" area. For the logio [ /([Nn]) ] 
Panel (c) of Figure 17 gives the spatial distribution of the 
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Figure 16. Panel (a): Diagnostic diagram of the logio [ /(^[g°]^) ] logic [ /(^Nnj) ] ''elation for points strictly associated to the bright, 
central structure and prior to the subtraction of the foreground/background material (see text). Definitions for blue crosses and red filled 
circles are the same as in Figure 12a. Panel (b): Diagnostic diagram of the logio [^ /(^[s"]^) j vs. logio j^ ^j^^jj j relation for points strictly 
associated to the bright, central structure and following the subtraction of the foreground/background material (see text). Each point 
kept its original symbol used in Panel (a). The hatched area roughly indicates the zone, in the parameter space, where already cataloged 
shock-dominated objects can be found. 



post-subtraction points using the same color code as in Pan- 
els (a) and (b). 

Combining the post-subtraction results for all three di- 
agrams allows us to spatially identify those sections on the 
plane of the sky that suggest shock excitation. Panel (d) of 
Figure 17 combines all results found in the first three pan- 
els. Red dots that can be found in all three post-subtraction 
diagnostic diagrams are again colored in red. These points 
targeted regions with high probability of shock excitation 
in central IC 1805. Those red dots that can only be found 
in Panel (c) are colored in green in Panel (d). These can 
be referred as interesting candidates for shock excitation al- 
though 2 out of 3 diagrams have failed to identify them. 
Blue dots in all three post-subtraction diagnostic diagrams 
are colored in blue and identify the zones, mostly located 
near the northern tip of the bright ionized feature, where 
photoionization clearly dominates. 

Shock-dominated material seems to prevail along the 
southern portion of the central structure, directly exposed 
to the stellar winds and currently eroded by the ionizing 
flux of the two most massive stars identified in Melotte 15 
(identified by black arrows in Figure lid). These are cat- 
aloged bv lMassev eFon (119951 ) as HB 15558 and HB 15570, 
two evolved 04111(f) (aaooo = 02''32™43', (52000 = 61°27'22") 



and 04If (qzooo = 02''32'"50', (52ooo = 61°22'42") stars re- 
spectively. Ionization by shocks clearly dominates the bright 
ionized gas at the periphery of the weak molecular material 
found in the southernmost portion of our FOV. Our results 
reveal that shock excitation appears to be correlated with 
bright ionized material, suggesting high electron densities if 
we assume the nebula to be optically-thin to the ionic light. 
This would then be easily correlated to the compressive na- 
ture of the hypothetical shocks. 

Extracting the electron densities p for all 1 419 points 
(rather in red or green in Figure lid) yields the blue-line 
histogram of Figure 10a. Its mean is estimated at 1251 50 
cm""', slightly below the mean value previously derived for 
the whole sample of 3 057 emission-line profiles retained for 
this work (see § 4.2.3). This hardly makes sense consider- 
ing the compression expected from shock waves. However, 
keeping in mind that the subtraction of, what we referred 
to as, the foreground/background material was required in 
order to spectrally detect evidences for shock excitation in 
IC 1805, it seems logical that the same procedure should 
be applied in the extraction of the electron densities specif- 
ically associated to the shocked, ionized component. Using 
the 1 419 post-subtraction emission-line profiles that showed 
evidences for shock excitation, the corresponding electron 
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Figure 17. Peak-intensity maps of the Hq!A6563 A ionic transition in IC 1805 (see caption of Figure 4). Panel (a): Obtained from 
the post-subtraction logio ^ /(^ [ g n ] ^) ] ^°§10 [ /(In n] ) ] relation, the specific lines-of-sight, associated to the bright, central structure and 
suggesting evidences for shock excitation, are represented in red. Blue dots indicate the position where the ionized material is most likely 
dominated by photoionization effects. Panel (b): Same as Panel (a) for the logio ^ i^ l s"]) ] [ s 'i j A6731 ratio relation. Panel (c): 

Same as Panel (a) for the logio [ i([n nj) ] }s !!] A6731 "'^'^ ratio relation. Panel (d): Combination of Panels (a) to (c). Red dots present 
in all first three panels are again colored in red. Dots colored in red in Panel (c) but blue in Panels (a) and (b) are colored in green. Blue 
dots present in all first three panels are again colored in blue. The two most massive O-stars, in the Melotte 15 cluster, are identified by 
the black arrows (see text). 



densities {= p shock) were computed following the method 
described in § 4.2.3. The red histogram, in Figure 10a, re- 
veals the distribution obtained for pshock- The mean elec- 
tron density along the line-of-sight at the position of shocks 
is estimated at 175 cm~^, a mean roughly 20 cm~^ greater 
than what was found for the black histogram in Figure 10a. 
This suggests, as expected, that the foreground/background 
material probably occupies, near Melotte 15, a large vol- 
ume filled with a very diffuse ionized component. This tenu- 
ous material undoubtedly contributes to attenuate the ionic 
emission of high-density condensations along the line-of- 
sight. 

Using the 1419 spectra of the foreground/background 
emission, we extracted to corresponding densities 
{= Pdiffv.se)- As expected, this tenuous component has 
a mean electron density, along the line-of-sight, of 25 
cm~"^. The distribution shows a narrow width of roughly 
15 cm~^ which reflects the homogeneity (in density) of 
the foreground/background material over the entire spatial 
extent of the central, ionized structure. Figure 18 shows 
the distribution of the pixel-to-pixel compression factor 
i.e., Pshock/ Pdif fuse- AH remaining 1419 points formed the 
black histogram. The mean compression is estimated at 8.5 



although the histogram is clearly dominated by values of 
pshock/ pdif fuse bclow 5. As cxpcctcd, the comprcssiou fac- 
tor is always greater than 1. Sub-distributions, in green and 
red, gives the compression factor respectively for the green 
and red dots obtained in Figure lid. The green histogram, 
formed of points suggesting possible shock excitation, is 
very similar statistically to the black distribution with a 
mean of 7.5. The red histogram, this time formed of points 
with a high probability of shock excitation, is slightly 
different, peaking at a greater pshock/ pdif fuse value. Its 
mean is estimated slightly above 10 i.e., those points (red 
in Figure 17rf) where, we believe, shocks may certainly 
dominate seem to be characterized by greater compression 
effects with respect to points (green in Figure 17ci) where 
the presence of shocks appears less obvious. 

For abiadatic jump conditions, Rankine-Hugoniot rela- 
tions indicate a compressi on factor of precis ely 4 between 
pre- and post-shocked gas (|Hennv et oZ.|[l999l . Chapter 12). 
The relatively large widths of the distributions displayed in 
Figure 18 indicate that isothermal shocks may dominate. 
Also, the fact that forbidden lines such as [N ii] and [S ii] are 
detected favors radiative cooling and therefore non-adiabatic 
effects [although these lines most likely originate from ion- 
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Figure 18. Histogram of the compression factor for points, associated to the central structure, revealing post-subtraction evidences for 
shock excitation (red and green dots in Figure ITd). The parameter pshock/ Pdif fv.se corresponds, pixel-to-pixel, to the ratio between 
the electron density computed from the [ g il j a6731 ratio following the subtraction of the foreground/background material and the 
electron density of the diffuse component retrieved from the foreground/background spectrum (see text). The green and red distributions 
are subsets of the black histogram and correspond respectively to the sets of green and red dots found in Figure lid. 



ized material trapped in the cooling, shocked outer shell and 
therefore may not be related to the wind-blown bubble's cur- 
rent state of expansion i.e., energy- or momentum-driven] . 
Theo retically, for an isothermal shock l|Dvson fc Williams! 
Il98d . Chapter 6), the compression factor corresponds to the 
square of the shock's Mach number, A/o (see the abscissa 
in Figure 18). Therefore, for a speed of sound of 10 km 
s~^ in the H"*" medium (in agreement with our choice of 
7400 K for the electron temperature in IC 1805; see § 4.2.3), 
shocks with typical velocities between 15 and 30 km s~^ (i.e., 
Mq =2—10) may have led to shock excitation near Melotte 
15. This is in agreement with typical expansion velocities 
found for wind-blown bubbles (see § 1). Shock velocities 
up to 50 km s~^ (i.e., A/o~25) are also found although 
such dynamics is certainly peculiar. Note, however, that the 
compression measured here could be an upper limit since 
photoeroded gas located in the molecular cloud's envelope, 
rather than the diffuse foreground/background component, 
may act as pre-shocked material. The cloud's envelope is 
expected to be, at least, a few tens of particles per cm"^ 
denser than typical va lues found for pdif fuse in this work 
(jRathborne et al\\2Q0^ . 

Assuming a coplanar geometry, roughly 3—4 pc sepa- 
rate the most massive star of the Melotte 15 cluster and the 
central, ionized structure. Usi ng the standard mod el for ex- 
panding wind-blown bubbles (|Weaver et fflZ.|[^977^ ■ the dis- 
tance reached by a given bubble's post-shocked shell can 
be estimated from the mechanical luminosity of the winds 
{Ln> = ■^Mu>v'^), the density of the pre-shocked medium 



{pdif fuse) and the total time during which winds have been 
blown by the central star {tu,)- Given tu, ~ 2.5 Myr, the age of 
the Melotte 15 cluster (see § 2), we expect that two massive 
giant and supergiant 04 stars have left just very recently 
the main-sequence branch. Therefore, for the calculations 
to follow, these were treated as 04V stars. From different 
studies proposing typical mass-loss rate s and terminal ve- 
locities with respect to spectral types (jPrinia et al ] |l990l : 
iLamers et a/.|[l999l : iMarkova et aZ.ll2004l ). and Voo were 
respectively estimated at 5 x 10~^ Mq yr~^ and 3 000 km 
s~^ for standard 04V stars. Given Pdi//iise~25 cm~^ (see 
above), nebular material found within 15—20 pc of Melotte 
15 is likely to have been previously disrupted by stellar winds 
emanating from the current cluster. This was calculated as- 
suming radiative (non-adiabatic) shocks ijLozinskaval Il993 . 
Equation 13.5). Given that interstellar shocks are usually 
adiabatic at the earliest times of the expansion, the non- 
radiative model yields a distance of roug hly 40 pc in the 
vicinity of Melotte 15 l|Lozinskaval Il992l Equation 13.9). 
Both distances exceed the extent of our FOV (15 x 9 pc^) 
and support the fact that shocks, attributed to the current 
star cluster, have had sufficient time to reach the bright 
structure in central IC 1805. 



5.2.2 Shock excitation in other zones of the nebula 

Besides the bright, central structure observed near the most 
massive stars of the Melotte 15 cluster, our FOV mostly 
reveals diffuse ionized gas, partially obscured by interstel- 
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lar dust. However, two structures are nonetheless detected, 
particularly well-defined in Ha (see § 4.1). In order to in- 
vestigate the presence of shock excitation in other portions 
of our FOV, a similar approach, to what was carried on in 
§ 5.2.1, is used here on these two features. 

5.2.2.1 South-east CO fragment. Panel (a) of Figure 
19 reveals the shape of the ionization front observed at the 
periphery of the CO fragment found in the south-east por- 
tion of our FOV (see Figure 8). As in Figures 3 to 9, North 
is up and East is left. The molecular feature is well-defined 
in the millimeter regime and its western side shows signs of 
erosion by the nearby star cluster. 

The ionized gas associated to the ionization front 
was first circumscribed while weaker areas near the bright 
H""" feature were selected and used to extract the fore- 
ground/background spectrum. In order to obtain a post- 
subtraction sample containing a sufficient number of 
emission-line profiles, only the first condition of § 4.2.2 was 
used. Therefore, results presented here should be cautiously 
interpreted. 

Prior to the subtraction, all points were initially found 
inside the "Hll regions" area of the logio ^ j^"]"*) j vs. 

lo 

Sio [ /(Inii] ) ] diagnostic diagram and black filled circles 
were used as symbols. Panel (b) displays the diagram 
for post-subtraction spectra only. These results reveal, as 
in § 5.2.1, that points are displaced toward the shock- 
dominated "SNRs" area although the tail is much less de- 
veloped when compared to Figure 166. No post-subtraction 
point actually enters the hatched area defined earlier. 

Shock excitation hence appears to have played a much 
more minor role in the ionization of this cloud's envelope 
when compared to the central structure of our FOV. The 
apparent shape of the ionization front indicates that the 
ionizing sources are located behind the molecular cloud 
(see § 4.1). Using kinematical information retrieved from 
the photoevaporated ma terial in the vicinity of the cloud, 
iLagrois fc Joncai (|2009al ) suggested that an appreciable dis- 
tance could exist between the CO fragment and the Melotte 
15 cluster (see the authors' Figure 11 and their associ- 
ated Flow G). Hence, large distances could signify here that 
stellar-wind shocks may have been partially dissipated be- 
fore they could reach the molecular cloud. Shock excitation 
is therefore likely measurable only within a certain distance 
of the shock sources i.e., stars with strong stellar winds in 
our case. Obviously, the stronger the shocks are, the greater 
this corresponding distance is. We reiterate that only the 
first 15 to 40 pc in the vicinity of Melotte 15 may have been 
disrupted by stellar winds (see § 5.2.1). At this point, the 
strongest wind shocks (fuelled up by the 04 stars) would 
have receded to the subsonic regime with velocities between 
4 and 9 km s~^ depending on the model used (radiative or 
non-radiative respectively). If, by any chance, the molecu- 
lar fragment found near the south-eastern boundary of our 
FOV is located relatively far from the ionizing sources (> 15 
pc), this could explain the apparent absence of strong, com- 
pressive shocks at its periphery. 

5.2.2.2 Cigar-like structure. Panel (c) of Figure 19 
shows the very tenuous, but well-defined, ionized coun- 
terpart of a cigar-like structure pointing toward a nearby 



star with strong Ha emission line. Again, North is up 
and East is left. The star is particularly visible on the 
left- hand side of Figu r e 8 and has been cataloged by 
iKohoutek fc Wehmevej (Il999i l as HBHa 6211-05 with a V- 
band magnitude of 14.6 and no known spectral type. Its 
coordinates are (aaooo = 02'"34"'10=.06, ($2000 = 61°24'35".7). 
The elongated feature was briefly introduced in § 4.1. 

Emission in Ha and [N 11] is detected on the outskirts 
of the cigar-like feature while its center appears mostly 
gas depleted. This overall scheme is particularly similar to 
what is expected fro m elephant trunks in Hll regions (e.g., 
ICarlavist et a/. !'2003Y Shadowing effects, caused by a dense 
neutral globule located in a radiation field, allow the warm 
gas behind it to recombine which could explain the absence 
of ionized material at the center of the elongated structure. 
The bright ionized rims are formed of photoevaporated ma- 
terial created by the erosion of the neutral globule. These 
fiows move away from the ionizing star creating the cigar- 
like shape. 

Both the FCRAO CO (1-0) survey and the Canadian 
Galactic Plane Survey at 21 cm (jNormandeau et aZ.lll997l ) 
reveal no indication for either molecular or atomic gas as- 
sociated to the elongated feature. Theoretically, H2 gas is 
expected to constitute the main component of the eroded 
globule while Hi material should result from recombina- 
tions in the tail. The whole structure, however, has rela- 
tively small angular dimensions (45" x 7") which suggests 
that emission from the neutral gas could be beam diluted in 
the low-resolution radio observations. 

Theory has suggested that shocks will deve l op in 
elephant trunks (|Mackav fc LimI I2OI0I : |Raga et al] I2OI0I ') 

although, in first approximation, the logio j^ ^^^^j^j j vs. 

logio [^ /(l^",]) ] diagram has shown nothing particular using 
the first set of Gaussian fits i.e., all points associated to the 
elongated structure are well-confined into the "H 11 regions" 
area. However, subtracting the foreground/background ma- 
terial yields the diagram displayed in Panel (d) of Figure 
19. None of the conditions listed in § 4.2.2 were considered. 
Only a S/N greater than 3 (sufficiently high to confirm a 
physical detection; see § 4.2.2) was here required. Therefore, 
the results are highly questionable although worth mention- 
ing. The ionized material at the center of the cigar-like fea- 
ture being very tenuous (or simply non-existent; see above), 
all points filling the diagram are found along the external 
bright rims. Panel (d) reveals that the "SNRs" area of shock- 
dominated material is highly favored, in agreement with the 
prediction of shock development proposed by the theory. 

The particular shape of the elongated structure pre- 
sented in Panel (c) could be easily mistaken for jets at- 
tr ibuted to Herbig-Haro (HH) objects. However, Figure 4a 
of lFrew fc Parkerl (|2010l ) indicates the peculiar low intensi- 
ties of the [N 11] lines in HH objects. This leads to a position, 

for HH objects in the logio [t(§^] vs. logio [i^l^] dia- 
gram, clearly above the "SNRs" area. This said, the diagram 
of Panel (d) most likely confirms that the cigar-like feature 
is not part of an HH object in IC 1805. Rather, the elephant- 
trunk nature, suggested here, appears very plausible. 
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Figure 19. Investigating tlie presence of shock excitation in other portions of IC 1805. Panel (a); Integrated intensity of the Ha ionic 
transition near the CO fragment, detected by the FCRAO survey, and located in the south-east portion of our FOV (see § 4.1). Panel 



(b): Diagnostic diagram of the logi 
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relation for points associated to the bright ionization front displayed 
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in Panel (a). The diagram is obtained following the subtraction of the foreground/background material (see text). Panel (c): Integrated 
intensity of the Ha ionic transition associated to a cigar-likc structure pointing toward a nearby star with strong H« emission. Panel (d): 
Diagnostic diagram of the logio [ /(^ [ sn/) ] ^°SlO [ i([n n ] ) ] ''^l^'tion for points associated to the cigar-like structure displayed in Panel 
(c). The diagram is obtained following the subtraction of the foreground/background material (see text). For Panels (a) and (c), both 
maps were processed using the spectral collapse of the initial mosaicked cube in the channel interval of the Ha A6563 A ionic line. 



6 CONCLUSION 

The use of the imaging Fourier transform spectrometer Spl- 
OMM allowed us to obtain series of emission-line profiles of 
the optical gas in the brightest, central portions of the Galac- 
tic IC 1805 H II region. The bandwidth used at data acquisi- 
tion allowed the simultaneous observations of the Ha A6563 
A, [N II] AA6548, 6584 A, and [S ii] AA6716, 6731 A ionic Unes 
(see § 3). 

The main goal of this work was to investigate on the 
presence of supersonic shock waves attributed to stellar 
winds in the vicinity of Melotte 15, the current star clus- 
ter actually fueling up the expansion of the IC 1805 nebula. 
Literature has long suggested that a kinematical detection of 
stellar winds, in H ii region, might represent a difficult task 
since the typical expansion velocity of shocked shells could 



be commonly confused with other dynamical processes re- 
vealing similar kinematical behaviors (see § 1). On the other 
hand, specific line ratios retrieved from the optical gas may 
indicate if the presence of ionized material is attributed to 
standard photoionization or to shock excitation. These line 
ratios therefore provide a non-kinematical tool for identify- 
ing shocks in the nebular volume. Our results are summa- 
rized as follow: 

I. The llj A(i'^!48 ^^^^'^ ratio genuinely deviates from the 
theoretical value of 3 (see § 5.1). Values varying between 
2.5 and 4 were commonly found. The distribution of the 
density measurements has not allowed to demonstrate that 
the [N ii] A6548 A line could be affected by collisional de- 
excitations, hence reducing its peak intensity. This scenario 
remains nonetheless plausible and could be verified if den- 
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sities, in the N"^ volume, could be measured precisely (see 
§ 5.1). Densities extracted from the [g "j a6731 ^^^^'^ may not 
perfectly reflect the conditions prevailing in the vol- 
ume if both nitrogen and sulfur are not co-spatial in central 
IC 1805. 

II. The d ia gnosti c diagrams, introduced by 
ISabbadin efai] (|l977l '). indicate, in first approxima- 
tion, that photoionization most likely dominate in IC 1805 
(see § 4.2.5). This initially holds even for the densest, most 
emissive structures in our field-of-view, directly exposed to 
the radiation fields and stellar winds of the nearby, most 
massive stars (see § 5.2). 

III. Evidences for shock excitation appear only follow- 
ing the subtraction of the diffuse foreground/background 
material. This component likely occupies a very large frac- 
tion of the nebular volume in IC 1805; a volume sufficiently 
large that, even though being tenuous in nature, the fore- 
ground/background material strongly dilutes the signal em- 
anating from shocked condensations along the line-of-sight. 
Its mean density is estimated at 25 cm~'^ (see § 5.2.1). 

IV. Oriented on a south-north axis and surrounded by 
numerous 0-stars, a bright, large ionized feature occupies 
the central area of our field-of-view. The last, tenuous frag- 
ments of an old molecular cloud can be found near its south- 
ern portion. Shocks may have contributed to ionize material 
found in the direct vicinity of the molecular clump while the 
northern parts of the ionized feature, deprived of molecular 
emission, appears to be largely dominated by photoioniza- 
tion (see § 5.2.1). 

V. The shock-excited ionized gas has a mean density of 
175 cm"'^ and the compression factor, between pre- and 
post-shocked gas, is typically between 2 and 10. For isother- 
mal shocks, this suggests shock velocities between 15 and 
30 km s~^, in agreement with models describing the ex- 
pansion of wind-blown bubbles (see § 5.2.1). Geometrically 
speaking, given the apparent proximity between the central, 
ionized structure and the most massive stars of Melotte 15, 
winds seem to have had sufficient time to reach the structure 
within a timescale corresponding to the age of the star clus- 
ter (see § 5.2.1). This gives credence to our assumption that 
shock excitation can be found in the south-central portions 
of our field-of-view. 

VI. Points identified with a high probability of shock ex- 
citation reveal compression factors typically greater than 
those points where shocks may be present although less cer- 
tain (see § 5.2.1). 

VII. Shocks did not seem to have played a major role in 
the ionization of a molecular cloud's envelope located in the 
south-eastern portion of our field-of-view. This most likely 
results from the molecular fragment being located too far 
from the ionizing sources so that (1) shocks induced by stel- 
lar winds have not reached yet the cloud or (2) shocks have 
reached the cloud with velocities too low to initiate shock 
excitation (see § 5.2.2.1). 

VIII. Shock development was clearly detected on the out- 
skirts of an apparently weak, but well-defined elephant 
trunk located in the eastern portion of our field-of-view (see 
§ 5.2.2.2). This is in agreement with theoretical works de- 
veloped on such feature typically found in Hll regions. 
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